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Abstract 



We extend the formalism developed in Ref. 20 for the renormalisation of 
the chargino-neutralino sector to the most general case of the MSSM with 
complex parameters. We show that products of imaginary parts arising from 
MSSM parameters and from absorptive parts of loop integrals can already 
contribute to predictions for physical observables at the one-loop level, and 
demonstrate that the consistent treatment of such contributions gives rise to 
non-trivial structure, either in the field renormalisation constants or the cor- 
rections associated with the external legs of the considered diagrams. We fur- 
thermore show that the phases of the parameters in the chargino-neutralino 
sector do not need to be renormalised at the one-loop level, and demonstrate 
that the appropriate choice for the mass parameters used as input for the 
on-shell conditions depends both on the process and the region of MSSM 
parameter space under consideration. As an application, we compute the 
complete one-loop results in the MSSM with complex parameters for the 
process ha — )■ xtx] (Higgs-propagator corrections have been incorporated 
up to the two-loop level), which may be of interest for SUSY Higgs searches 
at the LHC, and for chargino pair-production at an e^e~ Linear Collider, 
e+e" — )■ xtXj ■ We investigate the dependence of the theoretical predictions 
on the phases of the MSSM parameters, analysing in particular the numerical 
relevance of the absorptive parts of loop integrals. 
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1 Introduction 



The search for physics beyond the Standard Model (SM) is one of the main 
goals of the physics programme of the Large Hadron Collider (LHC). Su- 
persymmetry (SUSY) continues to be a particularly attractive extension of 
the SM. The minimal version, the Minimal Supersymmetric Standard Model 
(MSSM), predicts superpartners for all the fermions and gauge bosons of the 
SM as well as an extended Higgs sector consisting of two Higgs doublets. 

The recent signal discovered in the Higgs searches at ATLAS [Ij and 
CMS [2] with a mass of about 126 GeV, which is also compatible with the 
excess observed at the Tevatron |3|, is well in keeping with an interpretation 
in the MSSM in terms of the light (see e.g. [4]-[7]) or even the heavy CP-even 
Higgs boson fS^-^. On the other hand, the direct searches for superpartners 
at the LHC have so far not revealed any sign of a signal. The searches 
up to now have mainly been sensitive to the production of squarks of the 
first two generations and the gluino |9 , 10 . The analyses are just starting 



to become sensitive to the direct production of the squarks of the third 
generation and to the direct production of colour-neutral states of the MSSM, 



see e.g. Refs. 11 ,12 . In particular, in the chargino sector the most important 



limit is still the one from LEP of about 100 GeV 13 . The neutralino sector 
of the MSSM is even less constrained. If one drops the assumption of the 
GUT relation between the parameters M2, which appears in the chargino 
sector, and Mi, which appears only in the neutralino sector, the lightest 
neutralino can be arbitrarily light without violating the existing experimental 
constraints 
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In general, many of the MSSM parameters may take complex values. 
If this is the case, it leads to CP violation beyond that provided by the 
phase of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, i.e. the matrix 
that governs the mixing of the quarks in the SM. The CKM picture of CP 
violation has been remarkably successful in passing many experimental tests, 
which has led to stringent constraints on the possible structure of new physics 
contributions. On the other hand, additional sources of CP violation are in 
fact needed to explain the baryon asymmetry in the universe via a first order 
electro-weak phase transition, as the CKM phase can only partially account 
for the observed asymmetry (for a review see e.g. Ref. [Ts]). 

It is therefore of interest to investigate the collider phenomenology of 
CP-violating effects in the MSSM with complex parameters, see for instance 



Refs. 16 19 and references therein for studies of CP-violating effects at the 
LHC via rate asymmetries or triple products. While in some cases CP- 
violating effects occur already at tree-level, it is important to take into ac- 
count also loop-induced CP-violating effects. In order to make predictions 
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for CP-violating effects at the loop-level it is thus necessary to perform the 
renormalisation of the MSSM for the general case of complex parameters. It 
should be noted in this context that in the MSSM with complex parameters 
there are two sources of imaginary parts occurring at the loop level, namely 
complex pararameters and the absorptive parts of loop integrals. While ab- 
sorptive parts of loop integrals are often neglected in one-loop calculations, 
a consistent treatment of these contributions is essential in the MSSM with 
complex parameters since imaginary parts of the loop integrals can combine 
with imaginary parts of the MSSM parameters to contribute to the real part 
of the one-loop contribution. 

In this paper we work out the on-shell renormalisation of the chargino- 
neutralino sector of the MSSM for the most general case of complex param- 
eters, which involves in particular a consistent treatment of the absorptive 
parts of loop integrals. At leading order (LO) the chargino-neutralino sec- 
tor of the complex MSSM depends on the gaugino masses Mi and M2, the 
higgsino mass /i and tan /3, the ratio of the vacuum expectation values of the 
two neutral Higgs doublet fields, while at higher orders many more parame- 
ters become relevant. For the field renormalisation, we follow the formalism 



developed in Ref. 20 , where it was shown that it is convenient to choose 
different field renormalisation constants for incoming and outgoing charginos 
and neutralinos in order to ensure the correct on-shell properties of the exter- 
nal particles. While in Ref. 20 the parameters in the chargino and neutralino 
sector were assumed to be real, we perform the parameter renormalisation for 
the general case of complex parameters. We find that at the one-loop level 
the phases of the parameters in the chargino-neutralino sector do not need 
to be renormalised. Furthermore we demonstrate that an appropriate choice 
of the mass parameters used as input for the on-shell conditions depends 
both on the process under consideration and the region of MSSM parameter 
space. 



Recently, in Refs. 21 23 an on-shell scheme for the renormalisation of 
the chargino and neutralino sector in the complex MSSM has also been em- 
ployed. The differences in the approach of treating the CP-violating phases 
as compared to the present paper have been discussed in Ref. (24] (see also 



Refs. 25 26 



The on-shell renormalisation of the chargino-neutralino sector has been 



investigated for the MSSM with real parameters in Refs. 20 ,27 34 . More 



recently, in Ref. 35 a method to systematically choose on-shell conditions 
for the parameter renormalisation in any chosen scenario was presented, in- 
cluding the case of the parameters being strongly mixed. These results are 



in accordance with the earlier investigations of Refs. 20 , 34 



As an application of the framework for the renormalisation developed 
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in this work, we compute the complete one-loop results in the MSSM with 
complex parameters for the processes ha — )■ xtxj e"^e^ — )• xtxj- 
our numerical analysis we study in particular the dependence of the results 
on the phases of the complex parameters and we investigate the numerical 
impact of products of imaginary parts arising from complex pararameters 
and from absorptive parts of loop integrals. 

The decay of heavy Higgs bosons to charginos and neutralinos, ha — )• 
xtxj y is important in the context of SUSY Higgs searches at the LHC. The 
experimental signature of this process comprises four leptons and missing 
transverse energy |36j-|40j. In the search for the heavy neutral Higgs bosons 
of the MSSM this channel may provide sensitivity also in the "LHC wedge 
region" (see e.g. Refs. 41 , 42]) , where the standard searches for heavy MSSM 
Higgs bosons in t~^t~ and bb final states are expected to be not sufficiently 
significant for a discovery. Electroweak one-loop corrections to this class of 
processes have been evaluated in Refs. 43 45 for the case of real param- 
eters. As explained above, we have obtained the complete one-loop result 
for the decay of a heavy MSSM Higgs boson in the general case of complex 
parameters (which gives rise to a mixing of all three neutral Higgs bosons 
to form the mass eigenstates). We have incorporated into our result Higgs 
propagator corrections up to the two-loop level. 

Since charginos and neutralinos are expected to be among the lightest 
supersymmetric particles and, as mentioned above, their mass range is only 
weakly constrained from SUSY searches at the LHC so far, the direct pro- 
duction of these particles via e~^e~ annihilation is of particular interest for 
physics at a future e~^e~ Linear Collider (LC). We focus here on chargino 
pair-production, e+e~ — )■ xfxj- High-precision measurements of this pro- 
cess in the clean experimental environment of an e^e^ LC could be crucial 
for uncovering the fundamental parameters of this sector and for determining 
the nature of the underlying physics. Based on a leading-order treatment, 
the determination of the parameters Mi, M2, n and tan/3 is expected to be 
possible at the percent level, providing also sensitivity to non-zero phases 



of complex parameters 46 . At this level of accuracy, higher-order correc- 
tions need to be incorporated. For the case of real parameters, one-loop 
corrections to chargino pair production at a future LC have been investi- 
gated in Refs^j47jj48j . Loop-induced CP- violating effects have been studied 



in Refs. 49 , 50 , in particular the effect of complex parameters on certain 



asymmetries, but the considered quantities were UV finite and therefore did 
not require renormalisation. We extend the previous results to the general 
case of complex parameters. 

In the following section we will introduce the MSSM with complex pa- 
rameters and briefiy summarise our renormalisation procedure for the Higgs, 
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(s)fermion and gauge boson sectors. In Sec. [s] we present a framework for the 
on-shell renormalisation of the chargino and neutrahno sector of the MSSM 
for the general case of complex parameters, where in particular the imag- 
inary parts arising from complex parameters and from absorptive parts of 
loop integrals are consistently treated. As an application of this framework, 
in Sec. |4] we derive new results for two phenomenologically interesting pro- 
cesses, Higgs decays into charginos and chargino pair production at a future 
LC, and we study the dependence of the results on the complex parameters 
and the numerical impact of products of imaginary parts. We conclude in 
Seem 



2 The MSSM with complex parameters and 
its renormalisation in the on-shell scheme 

As mentioned above, complex parameters arise naturally in SUSY, inducing 
CP violation. In the most general MSSM (for the case of massless neu- 
trinos) there are 40 possible phases. Under the assumption, however, of 
Minimal Flavour Violation, which is motivated by the strong constraints 
on supersymmetric contributions to flavour-changing neutral current pro- 
cesses such as 6 — )■ S7, the number of MSSM parameters that may be com- 
plex reduces to 14: the phases of the sfermion trilinear couplings 0/, where 
/ = u,c,t,d, s,b,e, fi,T; the phases of the gaugino mass parameters (pMi, 
where i = 1,2, 3; the phase of the Higgsino mass parameter 0^, and the phase 
of the Higgsino mass mixing parameter 0^12 Out of these, the freedom 
to redefine fields means that any two phases may be rotated away, and as 
in Refs. 



20,52 



we choose those to be 0Af2 and 0rn.i2- We therefore consider 
12 non-vanishing phases in the following. The most restrictive experimental 
constraints on those phases arise from bounds on the electric dipole moments 



(EDMs) of the neutron dn, mercury dug and Thallium, dxi 53 55p" 

\dn\ < 2.910-2''ecm(90%C.L. 
I^Til < 9.010"25gpj^(^9Q%Q L_ 

I^Hgl < 3.110~2^ecm(95%C.L. 
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These EDMs are functions of CP-odd operators, e.g. the electron EDM d, 



60 



the quark EDM dq, the chromo EDM dq (for a recent review see Ref. 
via atomic or hadronic matrix elements which can have large theoretical un- 
certainties, ranging from ~ 10% to 50% 57,61 . Dominant contributions to 



^In addition, EDMs of the heavy quarks 56 , the electron 57p8 and the deuteron 59 
can have a significant impact. 
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dn, duo- and c/ti come from de, dud, dud- The MSSM contributions to these 



EDMs have been studied in detail in the hterature, see Refs. 60 ,61 and ref- 
erences therein. The dominant contributions to those EDMs involve the first 
two generations of squarks and sleptons, thereby imposing severe constraints 
on Aq^i for q = u,d,s,c and / = e, /x In our numerical evaluation below 
we set the severely constrained phases of the trilinear couplings to zero. In 
contrast, the third generation trilinear couplings of the squarks and sleptons 
are much less constrained by the EDM's, and can possibly result in large 
effects on observables. The phase 0^ of the higgsino mass parameter is also 
severely constrained in the convention where is rotated away. However, 
since it is the only phase present at tree level in the chargino sector, we will 
nevertheless investigate below the numerical impact of varying this phase. 
The bino phase (pMi , on the other hand, is less constrained by the EDMs, so 
that variations of this phase can potentially have interesting consequences 
in the neutralino sector. Cosmological effects of this phase include favouring 



bino-driven electroweak baryogenesis 64 , and modifying the relic density as 



well as both direct and indirect detection rates of neutralino dark matter. 
Concerning collider phenomenology, it can affect neutralino production rates 



and CP-violating observables at the LHC, see e.g. Ref. 61 , 65 , and at a 



future LC, see e.g. Refs. 16 52 66 68 



Below we will describe a systematic approach to the on-shell renormal- 
isation of the chargino and neutralino sector of the MSSM with complex 
parameters. In order to illustrate this approach, we will derive new one-loop 
results for two processes involving external charginos, namely ha — )■ xtxj 
and e~^e^ — )■ xtxJ, where i,j = 1,2 and a = 2,3. Besides the renormal- 
isation of the chargino and neutralino sectors, these processes also require 
the renormalisation of the MSSM Higgs sector, the sfermion sector, as well 
as the gauge boson and SM-fermion sector. We therefore first describe the 
renormalisation of these sectors before turning to the renormalisation of the 
chargino and neutralino sectors in the following section. 



2.1 Renormalisation in the Higgs sector 

The renormalisation of the Higgs sector plays an important role in the fol- 
lowing, firstly because we study the process ha — )■ xtxJ , where the external 
Higgs must be renormalised, and secondly because it enters the chargino- 
neutralino sector through the renormalisation of the parameter tan/3. Ex- 
pressing the Higgs potential in terms of the soft masses mi, m2, and the mass 



^Possible circumstances under which those bounds are evaded are discussed in Refs. 62 
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mixing parameter rriu leads to 



V, 



H 



mlHl^Hu + mlH;,H2i-e' 



(m 



li 



H2iH2i) + 



1 

2' 



:9"\HuH2^\\ 



(2) 



where g and g' are the U{1) and SU{2)l couphngs, and the two Higgs dou- 
blets can be expressed through 



•Hi 
•H2 



H12 

H21 
H22 



Vl 



V2 



^(02 + %) 



(3) 
(4) 



with the vacuum expectation values of the two Higgs doublets, Vi and V2, 
and tan/3 = V2/V1. The MSSM Higgs sector is CP-conserving at tree-level, 
i.e. the phase of the parameter rriu and the relative phase ^ between the two 
Higgs doublets can be rotated away and vanish upon the minimisation of the 
Higgs potential, respectively. Expanding Vh in terms of the neutral fields 



01/2 and xi/2 as well as the charged fields 
mass mixing terms, 

Vh = ... - T<^i0i - T0202 - T^iXi - ^X2X2 

/ 01 \ 

02 



results in tadpole terms and 



+ ^(01 02 Xl X2)M^^^^ 



Xi 

\X2 / 



1 



+ ... . 



(5) 



The two mass mixing matrices M,^,^^^ and M^±(^± can be diagonalised by 
rotation matrices parametrised by mixing angles a and /3„ and Pc, resulting 
in the neutral Higgs bosons, h, H and A, and the neutral Goldstone boson G, 
as well as the charged Higgs bosons if^ and the charged Goldstone bosons 
G^. By renormalising the Higgs doublet fields. 



Hi,2 ^ (1 + -6Zn,jn,,2, 



(6) 



one can obtain all the required Higgs field renormalisation constants, which 
can be written in terms of 6Zy^-^ ^, as discussed in Ref. 
to renormalise the tadpole coefficients. 
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. Similarly choosing 



\h,A — ^ Th,H,A + STh,H,A, 



(7) 
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the charged Higgs mass 

^M|± +5M|±, (8) 

and tan /3 via 

tan/3 tan/3(l + (5tan/3), (9) 

all parameter renormalisation constants can be obtained using relations con- 
necting them to STh,H,A, SM'^± and 6 tan (3, as also given in Ref. 69 . It 
is convenient to renormalise tan /3 in the DR scheme, see the discussion in 
Refs. 
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With the Higgs field renormalisation constants in the DR 



scheme, 



6Z. 



DR 

Hi 



-Re S 



/ div 

HH,a=0 



this yields 



5 tan (3 



DR 



DR 



(10) 

(11) 

(12) 



The mass of the charged Higgs is renormalised accoring to the usual on-shell 
condition, yielding 

(5M^± =Re Sh+h-(M^±), (13) 

and the renormalisation constants of the tadpole coefficients are fixed via the 
condition that the renormalised coefficient should vanish, leading to 

STh,H,A = -Th,H,A- (14) 

The numerically important Higgs propagator-type corrections in the MSSM 
Higgs sector not only affect the predictions for the Higgs boson masses, but 
also give rise to a loop-induced mixing between the neutral Higgs bosons. 
In order to ensure the correct on-shell properties of the external particles in 
the S-matrix elements, the mixing between different states has to vanish on- 
shell, and the residues of the propagators have to be normalised to one. We 
achieve this by applying finite wave function normalisation factors Zjj, which 
contain the complete one-loop contributions of the Higgs boson self-energies 
as well as the dominant two-loop corrections, as implemented in the program 



FeynHiggs 69 73 76 



The wave function normalisation factors Zjj, for which we use the defini- 



tion given in Refs. 79 ,80 , can be written as a non-unitary matrix Z. In this 
way a one-particle irreducible n-point vertex-function Gh^ involving a single 
external Higgs ha can be expressed as 



(15) 





(Gu 


J = z • 


Gh 


GhJ 


\Ga 
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in terms of the vertex functions Gh, Gh, Ga in the {h, H, A) basis. Here 
{ha,hb,hc) denotes some combination of {hi, h2, h^,). For definiteness we 
choose ha = hi, hf, = /i2 and h^ = h^. 

In addition to the mixing between the physical Higgs fields, a complete 
one-loop prediction for a process in the MSSM involving a neutral Higgs bo- 
son as external particle will in general also involve mixing contributions with 
the neutral Goldstone boson and with the Z boson. These contributions must 
exphcitly be included in the calculation at the one-loop level, as discussed in 



detail in Refs. 79 80 



2.2 Renormalisation in the sfermion sector 

As stated earlier, the calculation of e~^e~ — )■ xfxj at one-loop order requires 
the renormalisation of the sfermion sector, as enters the tree-level t-channel 
diagram. The renormalisation in the sfermion sector is furthermore needed 
for the evaluation of higher-order corrections in the Higgs sector, see above. 
At lowest order, the squarks and charged sleptons are mixed via 



Ml + mj + Mlili-Qfsl.) rufX} 
Mf=\ \ , 16 



for Sw as defined in the following subsection, making use of the abbreviation 
M| = M| cos 2/3, where Mz is the mass of the Z boson, and defining Xf by 

Xf = Af - fx* {cot (3, tan 13} , (17) 

where cot (3 applies for the up-type squarks, / = u,c,t, and tan/3 applies 
for the down-type sfermions, f = d,s, b, e, /i, r (we treat the neutrinos as 
being massless). Note that uif, Qf and // are the mass, charge and isospin 
projection of the fermion / respectively. Here M|, M? and the trilinear 
coupling Af are soft SUSY breaking parameters, of which only the latter 
may be complex, 

Af = \Af\e'''^f . (18) 

As we do not consider right-handed neutrinos, the sneutrino masses can be 
expressed by 

ml = Ml + ^Mlcos2(3, (19) 

where u = Ve, v^, Vt- 

In order to renormalise the sneutrino sector, we define the field and mass 
renormalisation constants by, 

v^{\^ ^6Z^)i) ml-^ml + 5m?. (20) 
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Imposing on-shell conditions in the sneutrino sector yields 



5Zf^ = -E'f^iml) and 5m? = Re E'^(m?). (21) 

It should be noted that choosing the sneutrino mass as an independent in- 
put parameter in this way implies that the renormalisation constant for the 
left-handed selectron mass is a derived quantity (following from SU(2) in- 
variance). 



2.3 Renormalisation in the gauge boson and fermion 
sector 

For the gauge-boson masses, Mw and Mz, we choose on-shell conditions. 
The weak mixing angle 9w is a derived quantity, following from 



2 



sin^^^^.^ = l-^. (22) 



With the renormalisation transformations 



Ml^Ml + 5Ml M^^M^ + 6M^, (23) 

sw ^ sw + Ssw, cw ^ cw + Scw, (24) 

where = cos^ 9w, this yields 

(5M^ = ReE^^(M^) and (5M| = ReE|^(M|), (25) 

where the transverse part Et(p^) of a self-energy T^^^{p) is defined according 
to 

S,.(P) = ( - + ^) St(/) - ^S,(/). (26) 
For the renormalisation of the weak mixing angle this results in 



^^w = 7, — -Jj2 ^72- = osw (27) 

2sw V J cw 

The renormalisation of the electric charge 

e ^ e(l + dZe) (28) 

in the on-shell scheme yields 



<5z, = ln-(o) + -^, (29) 



2 ' ' cw 



where ^^(0) 



ally 

constant, Aa, where 



. In order to avoid sensitivity to the light quark 
masses this is usually re-expressed in terms of the shift in the fine structure 



A a — n^ght fermions ( ^ ) ]\^2 ^ Tyight fermions ( )■ (30) 

z 

Here "light fermions" refers to the contributions of all quarks and leptons 
except the top quark. While the leptonic contribution to Aa can directly 
be calculated, the hadronic contribution is obtained from experimental data 
via a dispersion relation. In our numerical analysis below we will express our 
lowest-order results in terms of the fine structure constant at the scale M^, 
a(M|) =«/(! — Aa), so that the contribution of Aa is absorbed into the 
lowest-order coupling. 

Since we consider processes with external electrons, fermion field renor- 
malisation constants are needed. We define them according to the transfor- 
mation 

fL,R ^ (1 + \5Zl^)fL^R. (31) 

In the case where the mass of the fermion can be neglected, the on-shell 
condition leads to the simple expression 

5ZIj, = -S{^^(0), (32) 

where E{, are the left- and right-handed components of the fermion self- 
energy, respectively. 

While we do not consider processes with external gauge bosons, so that 
the field renormalisation constants of the gauge bosons drop out in our phys- 
ical results, these renormalisation constants do appear in expressions for 
individual vertices given below. For completeness, we therefore also list the 
expressions for the field renormalisation constants of the gauge bosons. We 
define the field renormalisation constants for the charged W bosons via 

^ {l + ]^5Zww)W^. (33) 
and for the neutral Z boson and photon by 

The on-shell conditions for the field renormalisation constants ensure that 
on-shell external particles have diagonal propagators with unity residues (see 
e.g. Ref. 81 and references therein). This leads to 

6Zvv = -(SD'(M2), (35) 
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where V = W, Z, 'y and S'(/c^) = ^^^2 ^ Ip2=fc2, and 

2S^^(M|) 2S^^(0) 
5^.z = ^^^7 = ^^^- (36) 

3 Renormalisation of the chargino and neu- 
tralino sector 

We now turn to the renormahsation of the chargino and neutrahno sector 
of the MSSM for the general case of arbitrary complex parameters. For the 
field renormalisation, we follow the formalism developed in the earlier work 



of Ref. 20 . We list here the relevant expressions for completeness. For 
the parameter renormalisation, we extend the results of Ref. |20j, which were 
restricted to the case of real parameters in the chargino and neutralino sector, 
to the general case of complex parameters. 

The charginos and neutralinos are the mass eigenstates of the gauginos 
and higgsinos, as seen from the relevant part of the MSSM Lagrangian, 

where cjl/r = 1/2(1 =p 75). The mass matrix for the charginos is given by 

^ _( M2 V2MwSi, \ 

^-^v^M^c, /i ^^^^ 

where sp/cp = sin/3/ cos /3. This matrix is diagonalised via the bi-unitary 
transform M + = U*XVl The mass matrix for the neutralinos in the 

{B, W, Hi, H2) basis is given by 

(Ml -Mzcpsw Mzspsw \ 



Y 



M2 MzCjjCw -MzSjjCw 

-MzCpSw MzC/sCw -fi 

\ MzSpSw -MzSpCw -f^ / 



(3^ 



where Mi is the bino mass. Since Y is complex symmetric, its diagonalisation 
requires only one unitary matrix A^, via M^-q = N*YN'^. The additional 
parameters that enter this sector are Mi, M2 and /i. 

In the MSSM with complex parameters, absorptive parts arising from 
loop integrals of unstable particles in general contribute to squared matrix 
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elements already at the one-loop level, since they can be multiplied by imagi- 



nary coefficients involving complex parameters. It has been shown in Ref. 20 
that a proper treatment of the absorptive parts from loop integrals of unstable 
particles implies that full on-shell conditions giving rise to vanishing mixing 
contributions on-shell can only be satisfied by the field renormalisation con- 
stants in the chargino and neutralino sector if they are allowed to differ for 
incoming and outgoing fields (see Ref. |88] for an earlier discussion of this 
issue in the context of the SM). Accordingly, we define the renormalisation of 
the chargino and neutralino fields in the most general way, i.e. we introduce 
separate renormalisation constants SZ^^^ and Z^^^ for incoming and out- 
going (left- and right-handed) charginos, respectively. The renormalisation 
constants for incoming and outgoing (left- and right-handed) neutralinos are 
denoted as ^Zq^^^ and Zq (^, respectively. Therefore, the renormalisation 
transformations for the chargino and neutralino fields read 

(^LXi ^ (1 + lSZ^)ijULXj, Xj^R -> X»~(l + \5Z^)ijUjR, 

^RXi ^ (1 + lSZ^)ijUjRXj, Xi^L -> X»"(l + \5Z^)ij^L, 

ujLxl ^ (1 + ¥Z^)ij^LX% ^iUJR ^ ^(1 + l5Z/^)ijUn, 

UJRX- ^ (1 + lSZ!^),,URx'j, ^ ^(1 + lSZ!^),,UL, (39) 

where the indices i,j can take values up to 2 for charginos and 4 for neu- 
tralinos, respectively. 

Concerning the parameter renormalisation, we treat Mi, M2 and fi as 
independent free parameters that are determined by imposing on-shell renor- 
malisation conditions in the chargino and neutralino sector. On the other 
hand the parameter tan /3 has been renormalised in the Higgs sector, and the 
parameters e, Mw and Mz (as well as the dependent parameter sin 9w) have 
been renormalised in the gauge sector, as described above, see Sec. [2j The 
renormalisation transformations for Mi, M2 and fi read 

|Mi I ^ |Mi I + 5|Mi 1 , 0Mi ^ 0A/1 + S(Pm, , 

IM2I ^|M2|+5|M2|, 

(f)^,-^ (f)f, + S(j)^, (40) 

where we treat the general case of complex parameters in the chargino and 
neutralino sector. As mentioned above, we have adopted the convention 
where the phase is rotated away. The parameter renormalisation in the 
chargino and neutralino sectors therefore amounts to the renormalisation 
of the five real parameters |Mi|, IM2I, (j)Mi and 0^. These parameter 
renormalisations induce a renormalisation of the mass matrices X and Y via 

X -^X + 6X, Y ^Y + 6Y. (41) 
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Our renormalisation scheme builds on and extends the work of Refs. 20 



|31]|34j|48]. Our scheme, besides addressing the general case of complex pa- 
rameters, differs from the methods followed in Refs. 



44,82 



where the renor- 
malisation was carried out for the case of real parameters. The approach 
of Ref. 82 differs from ours since in Ref. 82 the mixing matrices are left 



unrenormalised. In Ref. 44 the mixing matrices are renormalised using the 

Mz (and also 



83 , where the renormalisation of tan /3, Mi 



proposal of Ref. 

saiOy/) differs from our prescription. While the divergent parts of the pre- 



scription in Ref. [44] agree with the ones of the corresponding quantities in 

In an 



31 



31 



our approach, the finite parts differ, see the discussion in Ref. 
explicit comparison carried out for the case of real parameters in Ref. 
it was found that the resulting differences in the predictions for the physical 
chargino and neutralino masses based on the different methods were numer- 
ically small. 

As mentioned above, we determine the field renormalisation constants 
from full on-shell conditions that ensure vanishing mixing contributions on- 
shell for all chargino and neutralino fields. The propagators are required to 
have unity residues. The corresponding renormalisation conditions read 



2 = 
■m- "5 



lim — — 



O 9 



0, 



(42) 
(43) 



Xi 



where Xi = Xi ihJ = 1,^2) or = 1,2,3,4), and i ^ j. Here the 

renormalised propagator S^j^p"^) can be obtained from the IPI vertex func- 
tion f •^■'(p^), which can be expressed in terms of the renormalised self-energy 



Sf{p') = -(f If (p^))-\ f (/) = ^{^- m,)6,, + ^S,,(p2). 



(44) 



For convenience, we decompose the self-energy in terms of the coefficients 
Ej;/^(p2) and Sj.^/^^(p2) via 



(45) 



and define the left and right handed vector and scalar coefficients, S^-^^(p^) 
and Tj^^^^^{p'^) respectively, for the renormalised self-energy analogously. 



Note that the conditions in Eqs. (42) and (43) do not specify the wave- 



function renormalisation constants completely, and so in addition we impose 
the conditions that the renormalised propagators retain the same Lorentz 
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structure as the tree level propagators in the on-shell limit, i.e. 

t^ml) = tfMl), (46) 
«) = ti^ml). (47) 



Together Eqs. (42), (43), (46) and (47) result in the following expressions 
for the diagonal chargino wave function renormalisation constants, 

6Z^S = - ^S{ml.)-ml. [Si;,(m^.) + E£,(m^.; 

-Sg,(m^,) + iV5X^U'')u-{U*SXV%], (48) 
6Zi^^ = - T^Si^)-^ [Si;,(m^.) + S£,(m^.)]- 

-Sg,(m|±) + {V5X^U%,-{U*5XV%], (49) 

and the following off-diagonal chargino wave function renormalisation con- 
stants, 

Xi Xj 

-m^^ {U*5XV^)^^ - m,.^ {V5X^U^)^^] (50) 

Xj Xi 

, sr-SL/SR/ 2 N I ^SR/SL, 2 n 

-^4, - ^xf, iy^x^u^)^,] ■ (51) 

As a consequence of their Majorana nature, the renormalisation constants 
for the neutralinos satisfy the relations 

5<(f = ^Z^!,-- (52) 
The diagonal and off-diagonal wave function renormalisation constants are 
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^iven by 
5Z 



L/R 
OAi 



-S5:{fKo)-m|o[Sj;(m|„) + s£,(m|o)] 



(53) 



s:yL/R 



h|os5:{fKo)+m,om,oE-/;^Ko) 



2 2 



-m^o (N*6YN^)..-m^o (NSY^N^).]. 



(54) 



— r / D 

It should be noted that the barred constants, 8Z-J , are not related to 



(5Z 



via hermiticity relations, 



(55) 



For both the charginos and neutrahnos, SZ^J^ differ from {dZj^j^^Y in their 
absorptive parts only (arising from loop integrals of unstable particles), while 
this difference vanishes in the CP-conserving MSSM, see the discussion in 
Ref. 
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We briefly comment in this context on the similar issue arising in the 
fermion sector of the SM in connection to the renormalisation of the CKM 
matrix, see Refs. 



81,84-89 



It was realised that the hermiticity constraint 
for the incoming and outgoing fermion field renormalisation is incompatible 
with the demand of fulfilling the standard on-shell conditions. This is a con- 
sequence of absorptive parts of loop integrals, which are gauge-parameter 
dependent. Attempts to phrase the renormalisation prescription such that 
those absorptive parts do not enter the field renormalisation constants, restor- 
ing in this way the hermiticity relation, turned out to be problematic. The 
field renormalisation constants were found to be related via a Ward identity 
to the renormalisation constant for the CKM matrix, and the above prescrip- 
tion would lead to a gauge-dependent result 85 , 86 . Alternative methods 



to renormalise the CKM matrix have been proposed 85 , 87 , however in or- 



der to ensure that on-shell external propagators are flavour diagonal it was 
advocated to relax the hermiticity condition and to allow independent renor- 
malisation constants for incoming and outgoing fields 88 . This approach 

and 



was found to be consistent with the gauge invariance of the SM 
the CPT theorem 



89 



Due to the interference with the CKM phase, inclu- 
sion of the imaginary parts was found to give rise to numerically relatively 
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small shifts in the predictions for the relevant observables of ~ 0.5% 88 . 



In Sec. 4.1 we investigate the size of the effect of a proper treatment of the 
absorptive parts for the chargino and neutralino case and we analyse the 
numerical impact on predictions for physical observables. 

We now turn to the renormalisation of the parameters in the chargino and 
neutralino sector, which as discussed above comprises the renormalisation of 
the five parameters |Mi|, IM2I, |/i|, 0Mi and 0^. On-shell conditions for 
the parameters \Mi\, IM2I and can be obtained from the requirement 
that three of the chargino and neutralino masses are renormalised on-shell, 
in analogy to the case where Mi, M2 and fi are real, see Ref. [20]. For 
the two phases 0Mi and 0^ there exists no obvious on-shell condition (the 
same is true for several other MSSM parameters, for instance the parameter 
tan^). A possible choice would be to employ the DR scheme, as advocated 



in the "SPA conventions" 90 . However, we have verified explicitly that 
no renormalisation of 0Mi and 0^ is required at all in order to render the 
relevant Green's functions finite. This result can be understood as follows: 
starting from the Lagrangian expressed in the gaugino-higgsino basis, the 
diagonalisation of the mass matrices upon making the transition to the mass 
eigenstate basis leads to expressions in terms of the real combinations U*XV^ 
and N*YN''. In those expressions the phases of Mi and n that are present 
in the mass matrices X and Y have been compensated by the corresponding 
elements of the transformation matrices U, V and A^. Thus, the phases of 
Ml and n appearing in the couplings of neutralinos and charginos to other 
particles can be related to elements of the transformation matrices U, V and 
A^. The elements of those transformation matrices, however, do not need 
to be renormalised. This is in analogy, for instance, to the transformations 
of fields in the Higgs sector of the MSSM, where it is well known that the 
mixing angles a, (3-^ and (3^ (using the notation of Ref. [69]) do not require 
renormalisation. We therefore adopt a renormalisation scheme where the 
phases and 0^ of the parameters in the chargino and neutralino sector 
are left unrenormalised. This is convenient both from a technical and a 
conceptual point of view. 

We define the physical masses of the charginos and neutralinos Xi accord- 
ing to the real part of the complex pole 

^l = Ml-^M^.^x. (56) 

where F^. is the width of particle Xi- The physical mass at the one-loop level, 
M^., in general differs from the tree-level mass, m^^, by a finite amount. 
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Am^^., which is obtained from the relation 

M^, = m^,{l - ^Re [E,^(m|J + E,^(m|J]) - ^Re [Ef (m|J + (m|J] 
= m;^, + Am^,. (57) 

The renormahsation conditions for the three independent parameters |Mi|, 
IM2I and 1/^1 can be chosen such that three of the chargino and neutrahno 
masses are renormahsed on-shell, i.e. for those three particles the physical 
mass at the one-loop level is equal to the mass value at tree level, mj^^. 
Accordingly, these conditions can be written as 

Am^, ^ -^Re [Eij(m|J + E,^(m|J] - ^Re [Sf (m|J + E,f (m|J] 

= 0. (58) 

The resulting expressions for 5|Mi|, (5|M2|, depend on the choice that 
has been made for the three masses that are renormahsed on-shell. There 
are obviously three generic possibilities, namely selecting three neutralinos 
(NNN), two neutralinos and one chargino (NNC), or one neutralino and two 
charginos (NCC). Using the shorthands 

C(,) ^Re [m,+ [Ef:,,,(m^,) + E^,,(m^,)] +E^;<,(m^,) +E^5,(m^,)] 

- J] 25X,fcRe ([/,,• Vifc), (59) 



^ = 1,2 
fc = l,2 



iV(,) ^Re [m^o[So^,,,(m|,) + E«,(m|o)] + Sjt,(m|„) + Ej«(m|o)] 

-^45F,fcRe(%iV,,), (60) 



3=1,2 

fe=3,4 



the condition that the i'ih. neutralino mass is on-shell reads 

Ni, =2(5|M2|(Re (e-^<^^i A^^J + ReiV^J - 45|/i|Re {e-''>'>^Ni,^Ni,i), (61) 
while the condition that the i"th chargino mass is on-shell reads 

dn =2(5|M2|Re {Uin^V^.i) + 25\^l\Re {e-'^-Ui„2Vi''2). (62) 
In the NNN case, where neutralinos x°, and are chosen on-shell, we 



17 



obtain (5|Mi|, 5|M2|, 6\fi\ by solving Eq. (61) with i' = k simultaneously, 

Re {e-''''^NjsNji) Re Nl)Nk 



S\M,\ 



S\M2 



Q 



Q 



{Re{e-'^^NisNi4)ReN: 



+ (Re (e-*<^-iV,3iV,.4) Re Nj^ - Re {e-^'^-Nk^Nk,) Re N^N, 

+ (Re (e-*<^-iVfe3iVfc4) Re iV| - Re (e-^<^-iV,,3iV,4) Re Nl^)N^ , (63) 
1 

Q 

- Re {e-"*'' 
+ (Re {e-"* 



(Re {e-''''^Nj^Nji) Re (e'^^^^i A^fi) 



'iVfc3iVfc4)Re(e-*^"iiV2i) 



- Re (e- 
+ (Re (e 

- Re (e- 

1 



'iV.3iVi4)Re(e-''^*-^iiV2 



'iVfc3iV,4)Re(e-*<^«iiV/i))iV, 



2Q 



(Re A^l Re (e'^'^'^^i N]^) - Re (e"^"^"! A^"^!) Re N%)Nk 



+ (Re Nl Re (e'^^^i iV,\) - Re (e'^^^^i iV,\) Re iV22)iV, 
+ {Nl Re iV^^ - Re iV| Re (e-<^-i iV,\))iV; 
where 

: 2 (Re (e"*'^'^ A^is A^i4) Re A^^ Re (e-'^M, jys 



4\) 

Re A^l Re (e"*"^^ A^jsiVj^) Re (e"*"^*^! Nl^) 
Re (e"^"^^ Ar,3A/'i4) Re {e-'^^'^ Nf^) Re Ar22 



+ Re (e-*-^*^! A^fi) Re (e-^'^mj^Nj^) Re A^^ 
+ Re A^l Re (e-^-^^^iA^JJ Re (e-'^^'AT^gA^M) 
-Re (e-*<^*^iAr2 ) Re A^^^ Re (e-*^-Ar,3Ar,4)) . 



(64) 



(65) 



(66) 



For the NNC case, when neutralinos Xi^ Xj chargino xt on-shell, we 
solve Eq. (61) with i' = i,j and Eq. (62) with i" = k simultaneously, finding 



1 



5|Mi I = - ( 2 (Re (e-^'^'^ Ar,3Ari4) Re AT, 



i2 



- Re (e-*^«2iv|) Re (e-^'^''Ar,3Ar,4))Cfc 

+ (2Re (UMi) Re (e-*^-Ar,3Ar,4) 

+ Re ie-"^-Uk2Vk2)ReNf^)Ni - (Re (e 

+ 2Re (UkiVki) Re (e-^*-Ar,3Ar,4)) AT,- 



U,oV,.)ReN, 



i2 



(67) 
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5|M2| = ( - 2[Re Ni^Ni^) Re (e-*<^^^i AT 

- Re (e-*^-iV,3iV,4) Re (e"*^"! iV^ )]Cfc 

- Re (e-'^''t/fc2V^fc2) Re {e-'^^'^ Nf,)Ni 



2 1 



+ Re (e-^*-[/fc2Vfc2)Re (e-'^^iiV^ )iV,. j , (68) 

= ( - (Re iV^Re (e-^<^-i iVj,) - Re iVj^Re (e"^^"! iV^ )) 
+ Re (f/fciVfci) Re {e-"^"^ N^Ni 

- Re {UMi) Re (e-*<^-i iVfJiV,) (69) 
where 

R = 2Re (e-^^'^"f/fe2V^fc2) ( - Re iVf2Re (e"'^*^! Nf^) + N^Re N%) 
+ 4Re {UkiVki) ( - Re {e-^'^^NM^) Re (e-"^^^iiV2^) 
+ Re (e-'-^^iAT^) Re {e-'^^Nj^Nj^)). (70) 



Finally, if masses of one neutralino and two charginos, ^ xti 

shell, corresponding to the NCC case, we solve Eq. (61) with i' = i and 

Eq. (62) with i" = j, k simultaneously, resulting in 

Sm\ = - ^^^^ ((2Re (e-^^iV.3iV.,)Re {U,,V,,) 

+ ReNlRe {e-'^-U,2V,2))Cu + (Re [U.^Vj^] Re {e~"^-Uk2Vu2) 
- Re {e-'^-U,2V,2) Re {UuiVki))Ni 
-{ReNlRe{e-"^-Uu2Vk2) 

+ 2Re (e-'*-iV,3iV,4) Re {UuiVki))C^ , (71) 
6\M2\ = ^ (Re (e-"^''f/,2V^,2)C^fc - Re {e~''^-Uk2Vk2)Cj\, (72) 
= - ^{Re{U^rV,^)Ck-Re{UuiVki)C^, (73) 

where 

S = 2(Re (t/,iVfci)Re (e-^'^"f/,2V^,2) - Re (f/,iV,i) Re (e-''^-t/fc2li.2)) • (74) 

In order to apply the above renormalisation prescription to a certain 
process, it is necesary to investigate which of the possible choices of the three 
masses that are renormalised on-shell is in fact appropriate and results in a 
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well-behaved renormalisation scheme. It is usually convenient to impose on- 
shell conditions for the external particles of the process under consideration. 
However, some more care is necessary in order to ensure that the imposed 
conditions are indeed suitable for determining the parameters \Mi\, IM2I and 
\fx\. This issue was investigated in detail in Ref. for the case of the CPX 
benchmark scenario and a higgsino-like variant of the CPX scenario. 

We define here the CPX scenario such that the parameters take the values 
Ml = (5/3)(s2,/c2^.)M2, MsusY = 500 GeV, A,,i = 900 GeV, 0a/i = 0, 
0^ = 0, = 7i"/2, 4>Af3 = 7r/2, (pAfi2 ~ ^- gaugino-like case we 

use M2 = 200 GeV and fi = 2000 GeV, whereas for the higgsino-like case we 
choose fi = 200 GeV and M2 = 1000 GeV. In our numerical example below 
we furthermore use tan /3 = 5.5 and Mh± = 132.1 GeV. 





NNN 


NNC 


NCC 


NCCb 


NCCc 


NCCb* 


NCCc* 


5|Mi| 


-1.468 


-1.465 


-1.468 


2517 


-3685 


-365.4 


-4.671 


5\M2\ 


-9.265 


-9.265 


-9.410 


-9.410 


-9.410 


13.23 


13.23 


<5|/i| 


-18.48 


-18.98 


-18.98 


-18.98 


-18.98 


-5.333 


-5.333 


Am^o 

Al 











2517. 


-3681 


-5.809 


-0.522 


Am^o 

A2 








-0.1446 





0.3560 





-0.4806 


Am^o 





-0.5018 


-0.5016 


-0.8447 





-354.9 





Am^o 

A.4 


0.3238 


-0.1775 


-0.1775 


0.6851 


-1.439 


-0.1734 


-0.1548 


Am-± 


0.1446 


0.1445 

















Am-± 
X2 


0.5012 





















Table 1: Finite parts of parameter renormalisation constants and mass cor- 
rections in GeV for the CPX scenario (gaugino-like case). The last two 
columns, denoted with an asterisk, show the results for a higgsino-like sce- 
nario. The parameters for both scenarios are given in the text. 



In Tab. [T] we show the finite parts of the renormalisation constants 5|Mi|, 
6\M2\ and S\fi\ for the gaugino-like case of the CPX scenario, using five 
different choices of parameter renormalisation: 

• NNN with X2 X3 on-shell 

• NNC with X2 X2 on-shell 

• NCC with x?? Xi and X2 on-shell 

• NCCb with Xi ^md X2' on-shell 
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• NCCc with X3, xf and xt on-shell 

Also shown are the resulting one-loop corrections to those masses that are 
not renormalised on-shell. For two scenarios, NCCb and NCCc, also the 
results for the higgsino-like case of the CPX scenario are displayed (denoted 
as NCCb* and NCCc* in Tah.Q. 

For the gaugino-like case of the CPX scenario one can see from Tab. [T]that 
NNN, NNC and NCC are all suitable schemes, giving similar (and relatively 
small) values for the finite parts of the three renormalisation constants and 
modest corrections (in this example at the sub-GeV level) to the masses. On 
the other hand, the NCCb and NCCc prescriptions yield a huge value for the 
finite part of 5|Mi| and correspondingly an unphysically large correction to 
the mass ttl^o. This is due to the fact that this (gaugino-like) scenario has the 
hierarchy \Mi\ < IM2I <^ which implies that the parameters |Mi|, IM2I 
and broadly determine the values of the masses of x?, X2/X1 ^/Jxti 
respectively. Consequently, since the prescriptions NCCb and NCCc do not 
use the mass m^o as input, the parameter \Mi\ is only weakly constrained, 
yielding an unphysically large value for its counterterm and correspondingly 
an unphysically large correction to m^o . Thus, in order to avoid unphysically 
large contributions to 5|Mi|, the only bino-like particle in this scenario, x\) 
should be chosen on-shell. More generally, the renormalisation conditions 
must be chosen such that they provide sufficient sensitivity to each of the 
three underlying parameters that are renormalised, |Mi|, IM2I and 

On the other hand in the higgsino-like scenario, where -C |Mi| < IM2I, 
the parameters |Mi|, IM2I and form the dominant component of the 
masses of X35 x\/X2 a^'i Xi/2/Xi') respectively. Tab. jlj shows that in this 
case the results in scheme NCCc are well-behaved, i.e. it yields moderate 
contributions to the counterterms and the masses. This is a consequence of 
the fact that the bino-like Xs has been chosen to be renormalised on-shell. 
Scheme NCCb, on the other hand, where X2 instead of Xs is renormalised 
on-shell, shows unphysical bavaviour since the parameter |Mi| is only weakly 
constrained. 

The comparison of the gaugino-like and higgsino-like scenarios in Tab. [T] 
illustrates that an appropriate choice of renormalisation prescription depends 
not only on the process in question but also on the considered scenario of 
parameter values. It is therefore in general not possible to make a choice 
of the three masses that are renormalised on-shell in such a way that this 
prescription can safely be applied to all possible parameter configurations]^ 

•^One might wonder whether the problems related to selecting three out of six masses for 
an on-shell renormalisation could be avoided by using the DR scheme where the predictions 
for all the physical masses receive loop corrections. However, such a scheme will in general 
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Instead, it is necessary to adjust the renormalisation prescription such that 
at least one of the three masses that are chosen on-shell has a sizable bino 
component, at least one has a sizable wino component and at least one has 
a sizable higgsino component. Failing to fulfill this requirement will result 
in renormalisation constants being essentially unconstrained, and therefore 
taking large unphysical values. This issue was recently discussed for the case 
of real MSSM parameters in Ref. 35 , where it was similarly argued that one 
bino-, wino- and higgsino-like mass should be set to be on-shell|^ 



4 NLO predictions at the LHC and the LC 

We now wish to utilise the above renormalisation framework to make NLO 
predictions for the LHC and LC. This serves the purpose of illustrating the 
possible sensitivity of collider observables to the details of the renormalisation 
procedure, in particular the treatment of the imaginary parts, as well as 
highlighting the dependence of such observables on the CP violating phases 
introduced in the complex MSSM, which has so far not been studied in the 
on-shell scheme at one- loop. 

Specifically we calculate T{ha xtxj) a{e^e^ — )■ xtXj) at NLO, 
including full MSSM corrections and allowing Aj., fi, Mi and M3 to be 
complex. The diagrams are generated and the amplitudes calculated using 



FeynArts 91-95 , which however requires the counterterms for the relevant 
couplings as input. We calculated these counterterms by renormalising the 
fields and parameters as described in detail in Sec. |3j Explicit expressions 
for the necessary MSSM counterterms are given in the following subsec- 
tions. FormCalc f96'-(98] then was used to calculate the matrix elements 



and LoopTools |96| to perform the necessary loop integrals. The loop inte- 



grals are regularised via dimensional reduction 99 -101 , which ensures that 
SUSY is preserved, via the implementation described in Refs. 96,102 . We 
assume a unit CKM matrix. 

For both processes, we present our results for the scenario given in Tab.|2} 
In light of the current LHC results [9||10j , we take the masses of the first two 
generations of squarks and the gluino to be at 1 TeV. As the bounds on 



lead to a situation where the mass value inserted for a particle at an external line will 
be different from the mass value of the same particle if it appears as an internal line of 
a Feynman diagram. Such a mismatch is problematic in view of a consistent treatment 
of infrared divergent contributions associated with external particles that carry electric 
charge or colour. 



In addition, in Ref. 35 the case of large mixing was discussed in detail, and it was 
found that there the most stable results were obtained in an NNC scheme where the mass 
of a wino-like chargino is chosen on-shell. 
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Parameter 



|Mi| 
IM3I 

Ml,, 



Value 



100 GeV 
420 GeV 
1000 GeV 
1000 GeV 
400 GeV 
1300 GeV 
0.118 



Parameter 



Value 



M2 
Mh+ 
tan /3 

Ml, 

\M 

nit 



200 GeV 
800 GeV 
20 

500-800 GeV 
500 GeV 
1000 GeV 
173.2 



Table 2: Table of parameters used in our numerical analysis, where Ag and 
Ai denote the common trilinear couplings for the quarks and leptons, respec- 



tively. as{Mz) and mt are taken from Refs. 103 and 104 respectively. 



the third generation squark masses are much less constraining, we consider 
between 500 and 800 GeV. Here M^^ = (Mi)g, = (M/^),^ denotes 
the soft-SUSY breaking parameters as defined in Eq. (16) for the first gen- 
eration squarks, etc. In view of the fact that the LHC up to now places 
hardly any constraint on the charginos and neutralinos, we choose relatively 
low values for the mass parameters, |Mi| = 100 GeV, M2 = 200 GeV and 
= 420 GeV, adopting a CMSSM-like scenario in the chargino and neu- 
tralino sector. As we will be considering ha — )■ xtX2' where a = 2,3, we 
choose Mh± such as to ensure that the masses of /;.2,3 are above the threshold 
for these decay channels to be open. In view of the prospects for observing 
this decay via a signature comprising four leptons and missing transverse 
energy we choose relatively low slepton masses, i.e. Mj_^ ^ = 400 GeV and 
Mj^ = 500 GeV (where the value of the ratio M^^ 2IM1, has an impact on the 
relative amount of electrons and muons in the final state as compared to tau 
leptons), as well as relatively high tan /3, i.e. tan/3 = 20. Although further 
reducing the slepton masses and increasing tan (3 would enhance the signal, 
the EDM bounds would be tighter, as discussed below. Using the program 
FeynHiggs 69,73 76 and taking the current theoretical uncertainties from 
unknown higher-order corrections ito account, we have checked the predic- 
tions for the MSSM Higgs masses arising from the parameters in Tab. |2} 
Although the chosen parameters give predictions for the light Higgs above 



114 GeV, respecting the LEP limits 77,78 , they are not in keeping with the 



recent discovery of a scalar resonance at the LHC [1]|2]. As the purpose of 
this paper is not to study the detailed phenomenological consequences of the 
presented renormalization scheme, we do not discuss this issue further, but 
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note that a compatible light Higgs mass could be achieved, for example, by 
decreasing At to 1050 GeV. 

We study the effect of varying the phases c^^t, 4>At^ 'Pa^, 4>Mi, 4>M3 and 
(pn (using the convention that M2 is real). As discussed in Sec. g the EDM 
bounds on these phases can be quite restrictive, and we therefore evaluate the 



predictions for the EDMs explicitly using CPSuperH2 . 2 105 , incorporating 
further two-loop contributions using 2LEDM 60 . We find that for M^g = 



600 GeV the approximate bounds on the phases are (pAt ^ 7r/6, (pMi ^ vr/50 
and 0^ < 7r/1000. For Mg^ = 800 GeV the phase (pAt of the trilinear coupling 
in the stop sector is essentially unconstrained. The phases of M3, At, A^ are 
also found to be unconstrained. In obtaining these values, we took into 
account that while the prediction for is robust, the prediction for 
involves atomic matrix elements which are only known up to a factor 2 to 3. 
While the EDMs are mainly sensitive to the phases of the trilinear couplings 
of the first two generations, the relatively large value of tan /3 = 20 results 
in a non-trivial bound on the phase (\)At from the mercury EDM if M^.^ is 
sufficiently light. For tan/3 = 10, on the other hand, any value of would 
be allowed by the EDM constraints. Similarly, due to the choice of tan /3 in 
combination with the relatively low values of the slepton masses, the bound 
on is rather tight. For tan /3 = 10, on the other hand, the bound on 0Mi 
would be ~ 0A./j < tt/IO and, for example, upon additionally increasing M^,, 
for z = 1,2,3, to a common value of 600 GeV, would be unrestricted. 
We will discuss below the numerical impact of varying the different phases 
in our results within the context of the bounds on the phases arising from 
the EDM constraints. 



4.1 Heavy Higgs decays to charginos 

The recent discovery of a light Higgs-like state is not sufficient to distinguish 
the hypothesis of a SM Higgs boson from the hypothesis that the new state 
belongs to an extended Higgs sector. For instance, in the decoupling re- 
gion of the MSSM the new state could be interpreted as the lightest neutral 
MSSM Higgs boson. This state behaves SM-like in the decoupling region, 
while the heavy MSSM Higgs bosons decouple from the gauge bosons. In 
this region the class of processes involving heavy Higgs bosons decaying to 
pairs of neutralinos and charginos are of particular interest, as they could 
provide experimental evidence for an extended Higgs sector. Detection of 
these processes at the LHC could be possible in the final state with four 



leptons and missing transverse energy 36 ,38-40 . A study in the MSSM 
with real parameters came to the conclusion that with 300 fb~^ it may be 
possible at the LHC to detect heavy Higgs bosons H oi A with masses up 



24 



to ~ 800 GeV at the Sa level [37] . This could cover part of the "LHC wedge 
region" (see e.g. Refs. |4l|[42] ), where the standard searches for heavy MSSM 
Higgs bosons in r+r~ (or bb) final states will not be sufficiently sensitive to 
establish a signal. 

In the general case of the MSSM with complex parameters, where the 
three neutral Higgs bosons mix to form the mass eigenstates, we calculate 
the decay widths r{ha — > xtxj) for the two heavy MSSM Higgs bosons, i.e. 
a = 2, 3. Since in the parameter region of sufficiently high values of Mh± 
where these decays are open kinematically, the two states /12 and are 
nearly mass-degenerate, it will experimentally be very difficult to distinguish 
between these two state^in the signature with four leptons and missing Et- 



This fact is also apparent from the analysis of Ref. 37 , which was restricted 
to the case of real parameters, as the distribution of events arising from H 
and A in this analysis did not show considerable differences. 

The tree-level three-point vertex function for the interaction of charginos 
with neutral Higgs bosons hk = {h, H, A} is 

"(Born , , /^R \ , , I i ^l5fc:^/^L 

A minus sign appears between the uj^ and ujl terms for the CP-odd Higgs 
states, i.e. 5a:3 = 1 for /c = 3 and zero otherwise. The couplings, C^[^, are 
given by 

'^Shk = ^ijhk = ^25^^*^'"=' ^^^^ 



where 



Cijhk =i{o.kUj2Vii + bkUjiVi2), 

h ={-Ca,-Sa,CpJ. (77) 

Here Sa = sin a, Cq = cos a, etc., and the matrices U, V have been defined 
in Sec. |3j The diagrams for these decays at tree-level are shown in Fig. [T]for 
the example of the final state xtX2 ■ 

The tree- level decay width for the two-body decay ^ xtxj can there- 
fore be written as 

r--(/.. ^ x.^x") = ((|c^4J' + l^kl')-«'^^-^-) «. 



^ Complementarily, in Higgs decays to r leptons it may be possible to use certain 
asymmetries involving the r polarisation to analyse the CP properties of the decaying 



scalar particle 106 ■ 108 
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%1 

12 

Figure 1: Tree-level diagrams for the decay of neutral Higgs bosons /i, H and 
A to charginos and X2 ■ 



where 

K{x,y,z) = {{x-y-zf -Ayzf/^. (79) 

As explained in Sec. |2| we ensure the correct on-shell properties of the mixed 
neutral Higgs bosons by the use of finite wave function normalisation fac- 
tors Zij, which contain universal propagator-type contributions up to the 
two-loop level. With the aim to investigate the effect of the genuine vertex 
contributions for this process, we will compare our full one-loop result to 
an improved Born result which incorporates the (process-independent) nor- 
malisation factors Zij. Accordingly, we define the improved Born result by 
summing over the tree-level amplitudes for the three neutral Higgs bosons 
hk weighted by the appropriate Z factor, 

^Imp Born ^ QBorn ^^2 y (gg) 

k 

As mentioned above, by definition the Z factors do not include contributions 
due to mixing with the neutral Goldstone boson or the Z boson, and therefore 
the relevant one-loop contributions of this type must be explicitly included 
in the calculation. 

The one-loop corrections therefore involve vertex diagrams, examples of 
which are shown in Fig. |2| and the self-energy diagrams involving the Z 
and the Goldstone boson, examples of which are shown in Fig. |3j These are 
calculated following the procedure outlined earlier. In order to obtain UV 
finite results at the one-loop level, the three-point vertex function defined 



at tree-level in Eq. (75) must be renormalised, i.e. we need to calculate the 
diagrams shown in the bottom row of Fig. |4]for the vertex diagrams, and we 
also need to renormalise the self-energy corrections, i.e. we need to calculate 
the diagrams shown in the upper rows of Fig. |4} 



The counterterm for the three-point vertex function defined in Eq. (75) 
is of the form 

- + ^d-lY^-^5C^^,^, (81) 
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Figure 2: A selection of one-loop vertex diagrams for the decay of neutral 
Higgs bosons h, H and A to charginos Xi ^-iid X2 ■ 



and the coupling counterterms are given by 



1 



sw 

iR/L 



1=1 



RjL 



,R/L^ 



where, in analogy to Eqs. (76) and (77) 

nR _ r'L'' — ^ 



^3) 



Here 5Ze and 5sw are defined in Eqs. (29) and (27) respectively, and for 



the chargino field renormalisation constants, given in Eqs. (49) and (51) 



we have dropped the "±" as our tree-level diagrams do not contain any 
neutralinos. Note that the parameter renormalisation which enters these 
renormalisation constants is performed in the NCC scheme where the relevant 
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counterterms are defined in Eqs. (71 ) to (73). Tlie Higgs field renormalisation 
constants 6Zhh, SZhH, etc. appearing in Eq. (82) are linear combinations of 



the field renormalisation constants SZy^, SZy^ given in Eqs. (10) and (11) 
as specified in Ref. 
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In order to account for the mixing of the neutral Higgs bosons, the result 
for the amplitude is obtained by summing over the contributions of the three 
neutral Higgs bosons, multiplied with the corresponding wave function nor- 
malisation factors Zij, and by furthermore adding the mixing contributions 
of the Higgs boson with the Goldstone boson and the Z boson. 



^Imp. 1- 



Loop 



XfXj ha 



G,Z,se / 2 



Here 

butions in the MSSM, 



contains both the tree-level and the complete one-loop contri- 

denote the loop-corrected masses. 



and Ml, Ml 



while refers to the tree-level masses, see Refs. |79lp0] for further details. 

Since the vertex contributions involving virtual pnotons give rise to IR 
divergencies, we add to the contribution to the decay width arising from 



Eq. ( 84 ) the corresponding bremsstrahlung contribution from diagrams with 
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Figure 4: Counter-term diagrams for the decay of neutral Higgs bosons h, 
H and A to charginos xf and X2 ■ 



real photon emission, as shown in Fig. |5j In this way we obtain the complete 
1-loop result for the decay width, 

We have compared our result for the general case of complex parameters 
to the existing result that was restricted to the case of real parameters. 



given in Ref. 44 , and obtainable via the package HFOLD 45 . There the 



renormalisation prescriptions used for the chargino and neutralino mixing 



matrices, for tan /3 (where the renormalisation condition given in Ref. 109 



has been applied) and for the charge renormalisation differ from the ones 



used in the present work. The numerical evaluation in Ref. 44 was carried 
out in the SPS la benchmark scenario fllO|. We find agreement between our 
result and the result of Ref. |44j within the expected accuracy. 

For the numerical analysis of our results, we first consider the effects 
of the phases 0At, Mi and /i for the parameters as defined in Tab. |2] (as 
mentioned above, we use the convention where the phase of the parameter 
M2 is rotated away)]^ In Fig. |6] we show our results for r(/i2 — ^ X1X2) 
a function of the phases (pAt, Mi and fi. The full result corresponding to 



^ We also investigated the effects of the phases (pA^ , (pAt 7 't'Ah for the scenario described 
m Tab. [2] and found the maximum relative deviation from the decay width when the phases 
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Figure 5: Diagrams for the decay of neutral Higgs bosons h, H and A to 
charginos Xi ^"^^ X2 with real photon emission. 



Eq. (85) is compared with the improved Born result based on the amplitude 



given in Eq. (80 ). As mentioned above, the improved Born result incorporates 
higher-order contributions not only from the calculation of the mass of the 
decaying Higgs boson, but we have also dressed the lowest-order result with 
the universal wave function normalisation factors Zy. As one can see in 
Fig. |6} the dominant contribution to the dependence of the decay width on the 
three phases (pAt, <Pmi and 0^ is already present in the improved Born result. 
Despite the fact that (pAt only enters at loop level, its numerical impact on the 
decay width is found to be very significant, which is a consequence of the large 
(Yukawa enhanced) stop loop corrections in the Higgs sector. Compared to 
the decay width at (pAt = 0, the full MSSM one-loop corrected decay width 
is modified by up to 11% and 40% upon varying cpAt for = 600 GeV 
and 800 GeV, respectively. Note that here the larger impact of the phase at 
= 800 GeV is due to a threshold effect, as the mass of the decaying Higgs 



are zero to be far below the percent level, i.e. a maximum of 0.03%, 0.04% and 0.03%, 
respectively. 
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lies very close to the mass of the stops. The deviations from the improved 
Born result amount up to 3% and 1% for Mg.^ = 600 GeV and 800 GeV, 
respectively. 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



0A, "Aa 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



(piJn <pf,l 71 

Figure 6: Decay width in GeV for h2 — )■ X1X2 &s a function of the phases 
(j)At (upper), 0Mi (middle) and 0^ (lower). The sohd (dashed) lines show the 
1-loop corrected (improved Born) results for Mq^ = 600 (blue, left) and 800 
(red, right) GeV. 

The effects induced by the phase Mi, on the other hand, are relatively 
small. The shift compared to the decay width with (pMi = and the de- 
viations from the improved Born result are both at the percent level for 
= 600 GeV, and at the sub-percent level for M^^ = 800 GeV. As for 
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(f)At, the phase of Mi only arises in the expressions for the decay width at 
loop level, however the lower sensitivity of the decay width to this phase as 
compared to (pAt is expected as the pertinent loops are not Yukawa enhanced. 

Variation of the phase fi can in principle give rise to larger effects on the 
decay width of up to 8% and 10%, and deviations of 2% and 1% from the 
improved Born result, for Mg^ = 600 GeV and Mg^ = 800 GeV, respectively, 
as it appears in the tree level couplings. However, taking into account the 
tight constrained on this phase from the EDM bounds discussed earlier, the 
impact of the phase variation on the decay width is reduced to the sub-percent 
level. The kinks seen in the bottom right-hand plot of Fig. [6] showing the 
dependence of the decay width on (p^ for Mg.^ = 800 GeV arise due to the 
crossing of the masses of h2 and h-^ at these points. 

Accordingly, the phase having the most important impact on the decay 
width r(/i2 — )■ X1X2) is (pAf Due to the prospective difficulty in resolving 
the decays of h2 and experimentally, we have further investigated the 
dependence of the sum of the two decay widths, T{h2 — )• X1X2) + r(/i3 — )■ 
XiX2)j on (pAf As shown in Fig. [7| the marked dependence on (pAt is also 
present for the sum of the two decay widths, giving rise to shifts of up to 9% 
and 48% for Mg^ = 600 GeV and 800 GeV, respectively. 




Figure 7: Sum of the decay widths in GeV for /12 — )■ X1X2 X1X2 
as a function of the phase (pAf The solid (dashed) fines show the 1-loop 
corrected (improved Born) results for Mg^ = 600 (blue, left) and 800 (red, 
right) GeV. 

We now turn to the impact of the renormalisation procedure on our final 
result, i.e. we investigate the numerical relevance of the consistent treatment 
of the absorptive parts, which, as discussed in the previous section, affects 
in particular the field renormalisation prescription. In Fig. [8] we plot 

r p nFuU pimp. Bom 

4r = ^r-^^ , (86) 

p J^Imp. Born ' ^ ^ 
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where we compare the result including the absorptive parts with an approx- 
imation where the absorptive parts are neglected. For unpolarised charginos 
in the final state, the proper treatment of the absorptive parts affects the 
decay width by, at most, 0.4%. In Fig. [s] we display the decay width for 
polarised charginos in the final state, T{h2 — )■ XilX2r)- ^^^^ 
numerical impact of the proper treatment of the absorptive parts can amount 
up to a 3% effect in the decay width. On the other hand, as expected, the 
effect is seen to vanish for the case of real parameters, i.e. = 0, vr. The 
spikes seen in these plots arise due to the fact that at these values of (pAt 
the masses of the /12 and bosons cross. The spikes are seen to vanish for 
example on changing Mg.^ to 520 GeV, as the Higgs masses no longer cross 
for any value of 0^^, shown in the lower row of Fig. |8} While there may 
be a chance to determine the polarisation of charginos through the angular 
distribution of their decays products, a detailed study of the prospects at the 
LHC is yet to be undertaken. 



4.2 Chargino production at a future Linear Collider 

As a second example, we now investigate chargino production at a Linear 
Collider, a{e~^e~ — )■ xfxj). High-precision measurements of this process 
in the clean experimental environment of an e^e^ Linear Collider could be 
crucial for uncovering the fundamental parameters of this sector and for 
determining the nature of the underlying physics. A treatment addressing 
the most general case of complex parameters is mandatory in this context. 

At leading order, in the limit of massless electrons, the process a{e^e~ — ?■ 
X1X2) is described by the two diagrams shown in Fig. [o] (there is one addi- 
tional diagram for the xtXi and ^2^X2 final states). The transition matrix 



element can be written as 111 



e2 



Ma^{e+e xtXj ) = ijQai3 [v{e+)^f,ujau{e )] [u{Xj )Y^fiv{xt)] , (87) 

in terms of the helicity amplitudes Qap, where a refers to the chirality of the 
e~^e~ current, (3 to that of the xtxJ current, and summation over a and (3 
is imphed. 

Qll = Sij + DzGlC^+^- r, 

Qrl = 5ij + DzGrC^+^- ^ 



Qlr= 6i,+DzGL (c^+^-S + D^^(c^ + -VC- 

QrR = ^ij + DzGr 



R 



yR 
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Figure 8: Ratio of the 1-loop corrected decay width to the tree- level decay 



width, as defined in Eq. (86), for /i2 — > XilX-zr ^ function of the phase 
of At, (pAti showing the effect of absorptive parts of the self-energies in the 
chargino field renormalisation constants. The solid (dotted) line indicates the 
result with (without) taking the absorptive part into account. Parameters are 



as in Tab. [2| except 
for the lower plots. 
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500 GeV for the upper plots and Mq^ = 520 GeV 



The ZxtXj and eu^xt couplings are given by 

1 



C^+-_ = -\/,i, (89) 



and Gl, Gr, Dz and Dj;, are defined via 



o2 _ 1 1 

_^w 2 n — 
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Here Dz and Di, refer to the propagators of the Z boson and sneutrino, 
respectively, in terms of the Mandelstam variables s and t, and we can neglect 
the non-zero Z width for the considered energies. The tree-level cross section 
in the unpolarised case is then obtained by summing over the squared matrix 
elements, 

,tree 



K ' [S, 171-+, m- 



647r2s2 



JdnY,\Map\^. (91) 




Figure 9: Tree- level diagrams for the production of charginos xt X2 
the LC. 



The one-loop corrections involve self-energy, vertex and box diagrams, ex- 
amples of which are shown in Fig. 10 The diagrams are calculated following 
the procedure outhned earlier. 

In order to obtain finite results at one-loop, we need to renormalise the 
couplings defined at tree-level in Eq. (87), i.e. we need to calculate the dia- 
grams shown in Fig fill This involves renormalising the 'jxfxj, ZxfxJ and 
euext vertices as follows. 



CwSw 



n=l,2 ^ 



(92) 

where the analogous right-handed parts are obtained by the replacement 
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Xi Xi Xi 




h xi xi 



Xi 



Xi 



e ^0 Xi 



"Xn 



Xi 



X2 



X2 



Figure 10: Example one-loop self-energy (upper), vertex (middle) and box 
(lower) diagrams for the production of charginos and at the LC. 



L ^ R. Furthermore, 



•'-^^^i sw V V sw 2^ '^^ ^ ^ 

-{5Z,, + 5Z? )p2 - 5ml - -i^i^Z,, + 5Z? ) j , (93) 

for Ui = i>e, z>^, and the counterterm contributions of the coupling factors 
are given by 

SC^,..^0. (94) 

Note that again for brevity we drop the ± for the chargino field renormalisa- 
tion constants. Using this prescription to renormalise the vertices we obtain 
UV-finite results. 
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Figure 11: Counterterm diagrams for the production of charginos xt and X2 
at the LC. 

As the incoming and outgoing particles are charged, in order to obtain 
an infra-red finite result one must furthermore include soft photon radiation, 
which introduces the dependence on a cut-off. Using the phase-space slicing 
method the full phase space for the real photonic corrections can be divided 
into a soft, a hard coUinear and a hard non-coUinear (IR finite) region, 

^brems ^ ^soft(^^) ^ a^S,^(AE, AO) + a^^^U^E, Ae). (95) 

Here the singular soft and hard coUinear regions are defined by £■ < AE and 
9 < AO, respectively. Accordingly, the full cross section at next-to-leading 
order, including also the virtual contributions, is given by 

^fuU ^ ^tree ^ ^virt+soft^^^^ ^ (T^^^i^AE, AO) + a'^ZtcoA^E , AO). (96) 

Since in our analysis we are particularly interested in the relative size of the 
weak SUSY corrections, it is useful to consider a "reduced genuine SUSY 
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cross-section" , as defined by the SPA convention 90 , where the numerically 
large logarithmic terms of the QED-type corrections depending on AE and 
the terms proportional to = logs /ml are subtracted in a consistent and 
gauge-independent way. Accordingly, our numerical analysis below is done 



for the quantity (see Refs. 47 90]) 



a^' = - ( + log ^^(L. - 1 + A,)) , (98) 



TT \2 

where is given by the coefficient of the terms in the soft photon correction 
involving AE that arise from final state radiation (ISR) and from the inter- 
ference between initial and final state radiation (IFI). To be more explicit, 
we express the soft photon contribution as a sum of initial state radiation 
(ISR), FSR and IFI, 

= & + a + €1) ■ (99) 

Defining 5isr etc., in terms of the soft photon integrals lij [sl] , 

rISR _ T , T _ 9 7" 

"soft ~ ^PlPl ' -'P2P2 ^-'PlP2 5 

'^soft = ^ 2 (Jpjfc^ + /p2fe2 - Ipik2 ~ ^P2kl) , (100) 

the quantity can be obtained by taking the coefficient of Att log -^^ — — in 

AFSR , xlFl 
"soft ^ "soft- 

In the soft limit, the photonic contributions can be factorised into ana- 
lytically integrable expressions proportional to the tree-level cross-section for 
a{e~^e~ — )■ xtxj)- our calculation carried out with FormCalc the con- 
tributions from soft photon radiation have been incorporated using the soft 
photon factor as given explicitly in Ref. (sT]. 

Restricting our general result for complex MSSM parameters to the spe- 
cial case of vanishing phases, we have compared with the results given in 



Refs. 47, 481, which were evaluated in the SPSla' benchmark scenario. The 



renormalisation prescription in Ref. 47 differs from the one used in the 



present work in the renormalisation of the chargino and neutralino mixing 
matrices as well as of the electric charge and tan/3. Furthermore, a different 



choice has been made in Ref. 47 for the masses chosen as input in the se- 



lectron / sneutrino sector (as a consequence, the sneutrino mass in Ref. 47 
receives a shift at the one-loop level, while we have chosen an on-shell con- 
dition for the sneutrino mass). On the other hand, in the special case of 
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real parameters our renormalisation prescription is the same as the one used 



in Ref. 48 , with the exception of the renormahsation of tan/3. We find 



numerical agreement within the expected accuracy with the results given in 



Refs. 47,48 



We now turn to the investigation of our results for the case of complex 
MSSM parameters. In particular, we study the relative size of the one-loop 



corrections, 

5a ((T™^^'^ - a^'""^) 



a 



(101) 



as a function of 0^^, 0^ and 0Afi, for a -y/i = 800 GeV LC. In Fig. 12, the 
dependence on each of the phases is seen to be quahtatively the same for 
= 600 GeV and 800 GeV. In the case of 0At, the dependence is sizeable, 
due to the Yukawa enhancement for the stop loops, leading to effects of up 
to ~ 12% for = 600 GeV and up to ~ 6% for M^.^ = 800 GeV. For the 
EDM-allowed regions of 0a/i, which enters only at loop-level, and fi, which 
is highly constrained, the numerical impact of the phase variations is rather 
small, at most ~ 0.2% for (pMi- Thus, particularly for low values of Mg^, the 
phenomenologically most relevant effect arises from varying the phase (pAf 

In Fig. 13 we plot 6a /a as a function of 0At, comparing the results of 
including and ignoring the absorptive parts of the loop integrals in the field 
renormalisation. The left plot shows the unpolarised cross section for e~^e~ — )• 
xtX2 y while the right plot shows the result for specific polarisation states of 
the produced charginos, e~^e~ — )■ XilX2r- The impact of properly accounting 
for the absorptive parts of the loop integrals in the field renormalisation 



is clearly visible in Fig. [13] The resulting difference amounts to up to ~ 
2%, which could be phenomenologically relevant at linear collider precisions. 
Thus, a consistent inclusion of the absorptive parts of loop integrals is not 
only desirable from a conceptual point of view, but can also give rise to 
phenomenologically relevant effects in the MSSM with complex parameters. 
Another feature that can be seen in Fig. 13 is the fact that for (pAt = (in the 
present case where all other phases are set to zero) the results of including 
and ignoring the absorptive parts coincide. This is as expected, confirming 
that in the case of real parameters the absorptive parts can be neglected, see 
Sec.SH 

We find that the numerical impact of properly incorporating the absorp- 
tive parts increases with an increasing hierarchy between the parameters 



Ml, M2 and /i. This is illustrated in Fig. [T4| where the cross section for 
e^e^ — 7- xfx2 is shown as a function of fi. The difference between the re- 
sults including and neglecting the absorptive parts is seen to increase for 
increasing /i. The behaviour of the results around /i = 480 GeV is caused by 
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Figure 12: 6(t/(T for e+e — )■ X1X2 ^ function of the phases (pAt (upper), 
(f)Mi (middle) and 0^ (lower row) for Mg^ = 600 (left) and 800 (right) GeV. 
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Figure 13: 5a j a for e"'"e~ — )■ X1X2 (left) and e"'"e~ — )■ XilX-ir (right) as a 
function of the phase (pAt > ^oi Mg^ = 500 GeV, including (solid) and ignoring 
(dashed) the absorptive part of the loop integrals in the field renormalisation. 



a threshold effect due to the sneutrino mass lying in this region. 

5 Conclusions 

We have derived a renormalisation scheme for the chargino and neutralino 
sector of the MSSM that is suitable for the most general case of complex 
parameters. We have put particular emphasis on a consistent treatment of 
imaginary parts, which arise on the one hand from the complex parameters 
of the model and on the other hand from absorptive parts of loop integrals. 
We have demonstrated that products of imaginary parts can contribute to 
predictions for physical observables in the MSSM already at the one-loop 
level and therefore need to be taken into account in order to obtain complete 
one-loop results. 

Concerning the parameter renormalisation in the chargino and neutralino 
sector, we have shown that the phases of the parameters in the chargino and 
neutralino sector do not need to be renormalised at the one-loop level. We 
have therefore adopted a renormalisation scheme where only the absolute 
values of the parameters Mi, M2 and /i are subject to the renormalisation 
procedure. In order to perform an on-shell renormalisation for those param- 
eters one needs to choose three out of the six masses in the chargino and 
neutralino sector that are renormalised on-shell, while the predictions for the 
physical masses of the other three particles receive loop corrections. We have 
demonstrated, using the examples of gaugino-like and higgsino like scenarios 
with complex parameters, that the appropriate choice for the mass parame- 
ters used as input for the on-shell conditions depends both on the process and 
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Figure 14: 5a /a for e'^e~ — )■ X1X2 ^ function of /i, for Mg^ = 500 GeV 
and (pAt = 7r/4, including (solid) and ignoring (dashed) the absorptive part 
of the loop integrals in the field renormalisation. 



the region of MSSM parameter space under consideration. In order to avoid 
unphysically large contributions to the counterterms and the mass predic- 
tions one needs to choose for the on-shell renormalisation one bino-like, one 
wino-like and one higgsino-like particle. We have provided full expressions 
for the renormalisation constants of \Mi\, IM2I and \fi\ for the case where Mi 
and jU can be complex (i.e., we have adopted the convention where the phase 
of M2 has been rotated away) and for all possible combinations of charginos 
and neutralinos being chosen on-shell. 

For the field renormalisation, the consistent incorporation of absorptive 
parts gives rise to the fact that full on-shell conditions, which ensure that 
all mixing contributions of the involved fields vanish on-shell, can only be 
satisfied if independent field renormalisation constants are chosen for incom- 
ing and outgoing fields. If instead one works with a scheme where those 
renormalisation constants are related to each other by the usual hermitic- 
ity relations, non-trivial corrections associated with the external legs of the 
considered diagrams (i.e., finite wave function normalisation factors) need 
to be incorporated in order to obtain the correct on-shell properties of the 
incoming and outgoing particles. 

Within the described renormalisation framework we have derived com- 
plete one-loop results for the processes ha — )■ xfxj (supplemented by Higgs- 
propagator corrections up to the two-loop level) and e~^e~ — )■ xtxj iii the 
MSSM with complex parameters. For both processes we have investigated 
the dependence of the results on the phases of the complex parameters. In 
particular, we have analysed in this context the numerical relevance of the 
absorptive parts of loop integrals. 

Concerning our results for heavy Higgs decays to charginos, ha — )■ xtxJ ^ 
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which may be of interest for SUSY Higgs searches at the LHC, we find that 
the phase variations have a significant numerical impact on the prediction 
for the decay width. In particular, varying the phase (pA^, which is so far 
almost unconstrained by the EDMs, can lead to effects of up to 40% in 
the decay width. We find that the impact of the absorptive parts in the 
field renormalisation constants is most pronounced for the case of polarised 
charginos in the final state, for which the impact of a proper treatment of 
the absorptive parts can amount up to a 3% effect in the decay width. 

For chargino pair-production at an e'^e~ Linear Collider, e+e" — )■ xfxj, 
we find that the dependence of the cross-section on the phase (pAt is siz- 
able, yielding effects of up to 12% in our example. The impact of a proper 
treatment of the absorptive parts in the field renormalisation constants turns 
out to be numerically relevant in view of the prospective experimental accu- 
racy of measurements at a future Linear Collider. We find effects of 2-5% 
in our numerical example. Our results for the one-loop contributions to 
chargino pair-production at a Linear Collider for the general case of com- 
plex MSSM parameters may also be of interest for investigating the accuracy 
with which the parameters of the MSSM Lagrangian can be determined from 
high-precision measurements at a Linear Collider, since in this context the 
incorporation of higher-order effects in the theoretical predictions, which lead 
to a non-trivial dependence on a variety of MSSM parameters, is inevitable. ^ 
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